Search for: All records

Candidate Dark Galaxy-2 (CDG-2) is a potential dark galaxy consisting of four globular clusters (GCs) in the Perseus cluster, first identified in D. Li et al. through a sophisticated statistical method. The method searched for overdensities of GCs from a Hubble Space Telescope (HST) survey targeting Perseus. Using the same HST images and new imaging data from the Euclid survey, we report the detection of extremely faint but significant diffuse emission around the four GCs of CDG-2. We thus have exceptionally strong evidence that CDG-2 is a galaxy. This is the first galaxy detected purely through its GC population. Under the conservative assumption that the four GCs make up the entire GC population, preliminary analysis shows that CDG-2 has a total luminosity of L_V,gal = 6.2 ± 3.0 × 10^6 L_⊙ and a minimum GC luminosity of L_V,GC = 1.03 ± 0.2 × 10^6 L_⊙. Our results indicate that CDG-2 is one of the faintest galaxies having associated GCs, while at least ∼16.6% of its light is contained in its GC population. This ratio is likely to be much higher (∼33%) if CDG-2 has a canonical GC luminosity function (GCLF). In addition, if the previously observed GC-to-halo mass relations apply to CDG-2, it would have a minimum dark matter halo mass fraction of 99.94% to 99.98%. If it has a canonical GCLF, then the dark matter halo mass fraction is ≳99.99%. Therefore, CDG-2 may be the most GC dominated galaxy and potentially one of the most dark matter dominated galaxies ever discovered. 

We present MArk-dependently THinned POint Process (Mathpop), a novel method to infer the globular cluster (GC) counts in ultra-diffuse galaxies (UDGs) and low-surface brightness galaxies (LSBGs). Many known UDGs have a surprisingly high ratio of GC number to surface brightness. However, standard methods to infer GC counts in UDGs face various challenges, such as photometric measurement uncertainties, GC membership uncertainties, and assumptions about the GC luminosity functions (GCLFs). Mathpop tackles these challenges using the mark-dependent thinned point process, enabling joint inference of the spatial and magnitude distributions of GCs. In doing so, Mathpop allows us to infer and quantify the uncertainties in both GC counts and GCLFs with minimal assumptions. As a precursor to Mathpop, we also address the data uncertainties coming from the selection process of GC candidates: we obtain probabilistic GC candidates instead of the traditional binary classification based on the color–magnitude diagram. We apply Mathpop to 40 LSBGs in the Perseus cluster using GC catalogs from a Hubble Space Telescope imaging program. We then compare our results to those from an independent study using the standard method. We further calibrate and validate our approach through extensive simulations. Our approach reveals two LSBGs having GCLF turnover points much brighter than the canonical value with Bayes’ factor being ∼4.5 and ∼2.5, respectively. An additional crude maximum-likelihood estimation and simulation study show that their GCLF TO points are approximately 0.9 mag and 1.1 mag brighter than the canonical value, with p-values of ∼10^−8 and ∼10^−5, respectively. 

CDG-1 is a tight grouping of four likely globular clusters in the Perseus cluster, and a candidate dark galaxy with little or no diffuse light. Here we provide new constraints on the luminosity of any underlying stellar emission, using Hubble Space Telescope/UVIS F200LP imaging. No diffuse emission is detected, with a 2 σ upper limit of F200LP > 28.1 mag arcsec^−2 on the 5″ scale of CDG-1. This surface brightness limit corresponds to a 2 σ lower limit of >0.5 for the fraction of the total luminosity that is in the form of globular clusters. The most likely alternative, although improbable, is that CDG-1 is a chance grouping of four globular clusters in the halo of the Perseus galaxy IC 312. 

Recent observations of caustic-crossing galaxies at redshift 0.7 ≲ z ≲ 1 show a wealth of transient events. Most of them are believed to be microlensing events of highly magnified stars. Earlier work predicts such events should be common near the critical curves (CCs) of galaxy clusters (“near region”), but some are found relatively far away from these CCs (“far region”). We consider the possibility that substructure on milliarcsecond scales (few parsecs in the lens plane) is boosting the microlensing signal in the far region. We study the combined magnification from the macrolens, millilenses, and microlenses (“3M lensing”), when the macromodel magnification is relatively low (common in the far region). After considering realistic populations of millilenses and microlenses, we conclude that the enhanced microlensing rate around millilenses is not sufficient to explain the high fraction of observed events in the far region. Instead, we find that the shape of the luminosity function (LF) of the lensed stars combined with the amount of substructure in the lens plane determines the number of microlensing events found near and far from the CC. By measuringβ(the exponent of the adopted power law LF,dN/dL = ϕ(L)∝(1/L)^β), and the number density of microlensing events at each location, one can create a pseudoimage of the underlying distribution of mass on small scales. We identify two regimes: (i) positive-imaging regime whereβ > 2 and the number density of events is greater around substructures, and (ii) negative-imaging regime whereβ < 2 and the number density of microlensing events is reduced around substructures. This technique opens a new window to map the distribution of dark-matter substructure down to ∼10³ M_⊙. We study the particular case of seven microlensing events found in the Flashlights program in the Dragon arc (z = 0.725). A population of supergiant stars having a steep LF withβ = 2.55_−0.56^+0.72fits the distribution of these events in the far and near regions. We also find that the new microlensing events from JWST observations in this arc imply a surface mass density substructure of Σ_∗= 54M_⊙pc⁻², consistent with the expected population of stars from the intracluster medium. We identify a small region of high density of microlensing events, and interpret it as evidence of a possible invisible substructure, for which we derive a mass of ∼1.3 × 10⁸ M_⊙(within its Einstein radius) in the galaxy cluster. 

ABSTRACT We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the Amuse framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in Flash. We simulate a suite of 104 M⊙ clouds at 0.0683 pc resolution for ∼2 Myr after the onset of star formation, with virial parameters α = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial α having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over ∼0.01 Myr time-scales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses ≲1000 M⊙ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation. 

ABSTRACT The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled magnetohydrodynamics and direct N-body star cluster formation code torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the q ≥ 0.1 detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semimajor axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.